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Abstract 



<^ , In this paper we give a new derivation of the quark-antiquark potential in the 



Wilson loop context. This makes more explicit the approximations involved 

and enables an immediate extension to the three-quark case. In the qq case 

we find the same semirelativistic potential obtained in preceding papers but 
ft, 
D . 
c~| , for a question of ordering. In the 3q case we find a spin dependent potential 



identical to that already derived in the literature from the ad hoc and non 

correct assumption of scalar confinement. Furthermore we obtain the correct 

form of the spin independent potential up to the 1/m 2 order. 
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I. INTRODUCTION 

The aim of this paper is twofold. 

First we give a simplified derivation of the quark- ant iquark potential in the context of 
the so called Wilson loop approach [1] in which the basic assumptions, the conditions for 
the validity of a potential description and the relation with the flux tube model [2] can be 
better appreciated. 

Secondly we show how the procedure can be extended to the three-quark system [3] 
obtaining consistently not only the static part (stat) of the potential but also the spin 
dependent (sd) and the velocity dependent (vd) ones at the 1/m 2 order. 

For what concerns the qq potential the result is identical to that reported in [4,5] (see [6] 
for the spin dependent potential) except for a problem of ordering of minor phenomenological 
interest : 
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where 



Kgt = ~- + *r, (1-2) 

3 r 



8 \m{ mi J V 3 r 



—2 Si • (r x pi) -S 2 ■ (r x p 2 ) 

m\ W 2 



1 (Acts a 

2 \3 r 3 r 

+ -^— ^[S 2 • (r x Pl ) - Si • (r x p 2 )] + 
m 1 m 2 3 r 6 



14 f 1 

+ qM "3 

mim2 3 lr d 



3 

—(Si -r)(S 2 -r) - Si • S 2 



3 



>(r)Si • S 2 } , (1.3) 



vS = 7^— \t-( shk + r h r k )p h lP k 2 \ - 

2miin 2 lor J w 

2 1 1 

- Yl ^^{ a r P?t)w - t W r Pit • P2t}w • (1-4) 

~[ brrij J 6mim 2 

Obviously in Eqs. (1.2)-(1.4) r = z x — z 2 denotes the relative position of the quark and the 

antiquark and Pjt the transversal momentum of the particle j, p^ T = (5 hk — f h f k )p k - where 
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f = (r/r); the symbol { }w stands for the Weyl ordering prescription among momentum 
and position variables (see Sec. IV). Furthermore in comparison with [5] the terms in the 
zero point energy C have been omitted, since they should be reabsorbed in a redefinition of 
the masses in a full relativistic treatment. 
For the 3q potential the result is 
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Again 1^7 = Zj — z; denotes the relative position of the quark j with respect to the quark I 
(j, I = 1, 2, 3) and r, = z^ — zm the position of the quark j with respect to a common point 
M such that X)? =1 r, is minimum. As well known, if no angle in the triangle made by the 
quarks exceeds 120°, the three lines which connect the quarks with M meet at this point 
with equal angles of 120° like a Mercedes star (I type configuration, see Fig. 1A). If one of 
the angles is > 120°, then M coincides with the respective vertex and the potential becomes 
a two-body one (II type configuration, see Fig. IB). Furthermore V? \ is the Laplacian with 
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respect to the variable z 3 - and now v^. = (5 hk — fj"Tj)v k where r,- = (Tj/rj). The quantity 
z M in (1.8) is given by 

{N~ l y^_i (piT-/rjmA I type configuration 

J_1 V 3 J/ 3 3) (1.9) 

Pi/mi II type configuration : z M = z ; , 

iV being a matrix with elements N hk = Y^=i j-{^ hk ~ f^f k ). 

Finally Eq. (1.7) properly refers to the I configuration case. In general one should write 

Kd R = - E As, • V -)VJS x p; • (i-io) 
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In comparing (1.10) with (1.7) one should keep in mind that the partial derivatives in z^ 
of V^t vanish due to the definition of M. 

We observe that the short range part in Eqs. (1.6)-(1.8) is of a pure two body type: it 
is identical to the electromagnetic potential among three equal charged particles but for the 
colour group factor 2/3 and it is well known. Even the static confining potential in Eq. (1.6) 
is known [7,1,3]. The long range part of Eq. (1.7) coincides with the expression obtained by 
Ford [8] starting from the assumption of a purely scalar Salpeter potential of the form 

a (n + r 2 + r 3 ) ft&ft , (1.11) 

but at our knowledge it has not been obtained consistently in a Wilson loop context before. 
Eq. (1.8) is new. It should be stressed that (1.11) corresponds to the usual assumption of 
scalar confinement for the quark-antiquark system. As well known from this assumption 
Vj^at an d Kd 7 result identical to (1.2) and (1.3), but V^ turns out different from (1.4). 

The important point concerning Eqs. (1.1)-(1.8) and (1.10) is that they follow from 
rather reasonable assumptions on the behaviour of two well known QCD objects, W q q and 
Ws q , related to the appropriate (distorted) quark-antiquark and three-quark "Wilson loops" 
respectively. 

For the qq case the basic object is 

W q - g = ^(TYPexp(ig£dx»A^x)^ . (1.12) 



Here the integration loop T is assumed to be made by a world line I\ between an initial 
position y x at the time t[ and a final one Xx at the time t{ for the quark (t^ < if), a similar 
world line 1^ described in the reverse direction from X2 at the time U to y2 at the time 
t; for the antiquark and two straight lines at fixed times which connect Xi to X2, y2 to yi 
and close the contour (Fig. 2). As usual A^(x) = |A a A"(a;), P prescribes the ordering of 
the color matrices (from right to left) according to the direction fixed on the loop and the 
angular brackets denote the functional integration on the gauge fields. 

The quantity i In W q q is written as the sum of a short range contribution (SR) and of a 
long range one (LR): i In W q q = i In W^ + i In W^. Then it is assumed that the first term 
is given by the ordinary perturbation theory, that is at the lowest order 

i In Wjf = \g 2 f dx'i f dx^D^( Xl - x 2 ) (1.13) 

[Ppuv being the usual gluon propagator and a s = g 2 /4:7r the strong interaction constant) and 
the second term by the so called "area law" [9,1,4] 

tlnW q f = aS min , (1.14) 

where S m i n denotes the minimal surface enclosed by the loop (a is the string tension). 
Obviously Eq. (1.13) is justified by asymptotic freedom, Eq. (1.14) is suggested by lattice 
theory, numerical simulation, string models and other types of arguments. 

Up to the 1/m 2 order, the minimal surface can be identified with the surface spanned 
by the straight line joining (t, Zi(t)) to (t, z 2 (t)) with U < t < tf, the generic point of this 
surface is [4] 

< in = t u min = sz 1 (t) + (1 - s)z 2 (t) (1.15) 

with < s < 1 and z x (t) and z 2 (t) being the positions of the quark and the antiquark at 
the time t. 

We further perform the so called instantaneous approximation in (1.13), consisting in 
replacing 
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-oo 

and use (1.15) and (1.16) at an early stage in the derivation procedure. In this way we shall 
obtain Eqs. (1.1)-(1.4) in a much more direct way and without the need of assuming a priori 
the existence of a potential as done in [4]. So, once that Eqs. (1.13) and (1.14) have been 
written, Eqs. (1.15) and (1.16) give the conditions under which a description in terms of a 
potential actually holds. 

Notice that, while (1.12), (1.14), and even (1.13) in the limit of large U — t- u are gauge 
invariant quantities, the error introduced by (1.16) is strongly gauge dependent. The best 
choice of the gauge at the lowest order in perturbation theory is the Coulomb gauge for 
which the above error is minimum. To this choice Eq. (1.4) does refer. 

For the three-quark case the quantity analogous to (1.12) is 
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Here dj,bj are colour indices; Tj denotes a curve made by: a world line Tj for the quark 
j between the times t\ and if (£; < if), a straight line on the surface t = t[ merging from 
a point / (whose coordinate we denote by dm) and connected to the world line, another 
straight line on the surface t = if connecting the world line to a point F with coordinate Xm 
(Fig. 3). The positions of the two points / and F are determined by the same rules which 
determine the point M above. 

The assumption corresponding to (1.13), (1.14) is then 

2 

r 

i<j 

where now S m i n denotes the minimum among the surfaces made by three sheets having the 
curves Ti, T2 and T3 as contours and joining on a line Tm connecting / with F. 

We shall see that Eqs. (1.5)-(1.8) follow if we substitute (1.16) in (1.18) and again replace 
Smin with the surface spanned by the straight lines 

6 



dn^g = o9 2 Y. dx i / drfiD^Xi -Xj) + aS min , (1.18) 



«? min = t llf = SZj (t) + (1 - S )z M (t) (1.19) 

with j = 1,2,3, s G [0,1], z M (t) being again the point for which Y%=i \ z j(t) ~~ z m(^)| is 
minimum. 

The plan of the paper is the following one. In Sec. II we shall report the simplified 
derivation of the quark-antiquark potential as sketched above. In Sec. Ill we shall report 
the derivation of the three-quark potential. In Sec. IV we shall make some remarks and 
discuss the connection with the flux tube model. 

II. QUARK-ANTIQUARK POTENTIAL 

As usual the starting point is the gauge invariant quark-antiquark (gi, q 2 ) Green function 
(for the moment we assume the quark and the antiquark to have different flavours) 

1 — — c 

G(x 1 ,X2,yi,y 2 ) = -(0\T^2( y x 2 )U( y x 2 ,x 1 )ij 1 ( y x 1 )ip 1 (y 1 )U(y 1 ,y 2 )4> 2 (y2)\0) = 

= ^Tr(U(x 2 ,x 1 )Sf(x 1 ,y 1 \A)U(y 1 ,y 2 )C~ 1 S!(y 2 ,x 2 \A)C) . (2.1) 

Here c denotes the charge-conjugate fields, C is the charge-conjugation matrix, U the path- 
ordered gauge string 

U(b,a) =Pexp lig f dx^ A^x)\ (2.2) 

(the integration path being the straight line joining a to b), Sf and S 2 the quark propagators 
in an external gauge field A^; furthermore in principle the angular brackets should be defined 

as 

fV[A]M f (A)f[A]e iS M 

m*\>- f V [ A ] Mf ( A ) e iS[A] > v- 6 ) 

S[A] being the pure gauge field action and Mf(A) the determinant resulting from the explicit 
integration on the fermionic fields. In practice assuming (1.14) corresponds to take Mf(A) = 
1 (quenched approximation). 



Summarizing the first part of the procedure followed in Ref. [4] (see such paper for 
details) first we assume x® = x^ = if, y\ = y® = U with r = if — 1 ; > and note that Sj are 
4x4 Dirac indices matrices type. Then performing a Foldy-Wouthuysen transformation on 
G we can replace Sj with a Pauli propagator Kj (a 2 x 2 matrix in the spin indices) and 
obtain a two-particle Pauli-type Green function K. We shall show that in the described 
approximations this function satisfies a Schrodinger-like equation with the potential (1.1)- 
(1.4). 

One finds (see [4]) that, up to the 1/m 2 order, Kj satisfies the following equation 

. d 



dx° 



K j (x,y\A) = H FVf K j (x,y\A):-- 



1 I A\2 * , A\4 9 

■oj -o-r 

Kj(x,y\A) (2.4) 



mj + — (P, - 9 Af - -^( Pj - gAf - J-Sj -B + gA 
Zrrij Srrij rrij 



^{diE 1 - ig[A\ &]) + ^e ihk Sf{( Pj - gA)\ E h ] 



with the Cauchy condition 

K j (x,y\A)\ 3 » =1 p = 8\*.-y) (2.5) 

where e lhk is the three-dimensional Ricci symbol and the summation over repeated indices 
is understood. By standard techniques the solution of Eq. (2.4), with the initial condition 
(2.5), can be expressed as a path integral in phase space 

Kj(x, y\A) = / V[zj, pj] Texp li dt [ Pj ■ ±j - H FW ] \ ; (2.6) 

J Zj (y°)=y { Jy° ) 

here the time-ordering prescription T acts both on spin and gauge matrices, the trajectory 
of the quark j in configuration space is denoted by Zj = Zj(t), the trajectory in momentum 
space by pj = Pj(t) and the spin by Sj. Then, by performing the translation 

p^p + gA, (2.7) 

we obtain an equation containing the expression dt (gA° — gz ■ A) = g dx^A^, which is 



formally covariant. [] It is also useful to have an expression for Kj in which the tensor field 
F^ u and its dual F^ v appear. To this end we make the further translation 
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and, apart from higher order terms, we obtain 
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Thus we obtain the two-particle Pauli-type propagator K in the form of a path integral 
on the world lines of the two quarks 
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K(xi,x 2 ,yi,y 2 ;r) = / P[zi,pi] / £>[z 2 ,p 2 ] 
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Here T s is the time-ordering prescription for spin matrices, P is the path-ordering prescrip- 
tion for gauge matrices along the loop Y and as usual 



F^ = d»A v - d v A* + ig[A*, A v \ 



(2.11) 
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(2.12) 



1 More precisely, since the A h are matrices, the step / d 3 p f(p — gA) = J d 3 p f(p) can be justified 
by expanding f(p—gA.) in powers of g; apart from the zeroth order term, all the other terms involve 
derivatives of /(p) and do not contribute to the integral. 
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D v F vli = ff , F vli + ig[A v 1 F l 
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(2.13) 



and e^ vpu is the four-dimensional Ricci symbol. 

Furthermore as in Eq. (1.12) I\ denotes the path going from (tj,yi) to (£f,xi) along 
the quark trajectory (£,zi(t)), T 2 the path going from (£f,x 2 ) to (£j,y 2 ) along the antiquark 
trajectory (£,z 2 (£)) and T is the path made by r\ and T 2 closed by the two straight lines 
joining (£;,y 2 ) with (£;,yi) and (tf,xi) with (tf,x 2 ) (see Fig. 2). Finally Tr denotes the trace 
on the gauge matrices. Note that the right-hand side of (2.10) is manifestly gauge invariant. 

What we have to show is that the angular bracket term in Eq. (2.10) can be expressed 
as the exponential of an integral function of the position, momentum and spin alone taken 
at the same time t: 
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V qq playing the role of a two-particle potential. To this aim, expanding the logarithm of the 
left-hand side of (2.14) up to 1/m 2 order, we should have 
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with the notation 



«/H» 



\(YY¥[eMig§rdx» A^x))}f[A}) 
\ (Tr P exp(ig § r dx^ A^x))) 



(2.17) 



and W q q defined in Eq. (1.12). 

At this point in Ref. [4] we assumed that a quantity V qq satisfying (2.16) existed and 
derived its form. Here we no longer make such an a priori assumption but start directly 
from (1.13) and (1.14). 
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1. Contribution to the potential coming from i\uW c 
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In the Coulomb gauge we have 
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where D F (x) = J -0- )k ke ikx . 

The pure temporal part Dqq{x) is already of instantaneous type, for the pure spatial part 
we have the instantaneous limit 
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Replacing (2. 18), (2.20) and (2.21) in (1.13) we have 
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Similarly, if we denote by u^ = M M (s,t) the equation of any surface with contour F 
(s G [0, 1], i G [ti, if], u°(s,t) = t, u(l,i) = Zi(i), u(0, i) = z 2 (i) ), we can write 
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where the index T denotes the transverse part with respect to the unit vector s: 
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Then in approximation (1.15) we have 
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<9u min 

-^— = Zl (t) - Zi(t) = t(t) , (2.25) 

<9u min 

= szi(t) + (1 - s)z 2 (t) (2.26) 



and so 
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(2.27) 



where obviously now s — f and we have expanded the square root and performed the s 
integration explicitly. 

Notice now that, at the lowest order, ij can be replaced by Pj/rrij in (2.22) and (2.27). 
Then such equations become of the correct form required by (2.16) and so does the entire 
i In Wgq. In conclusion we have a first contribution to V g7i (pure Wilson loop contribution) 
in the form V^it + Vvl with V^L and V™ as given in (1.2), (1.4). The ordering prescription 
in (1.4) shall be discussed in Sec. IV. 

2. Spin related potential 

To obtain the remaining part of the potential we must evaluate the expectation values 
of the form (2.17) occurring in (2.16). 

Let us consider an arbitrary infinitesimal variation zi(£) — ► z i{t) + Szi(t) vanishing at 
t — ti and t — £; and evaluate 6(i\nW q q). | From (1.12) we have 



2 The result is easily achieved in the case of an abelian gauge theory: in fact one can freely commute 
the fields A„. In our case, on the contrary, the gauge theory is non abelian: one can still commute 
fields referring to different points because of the presence of the path-ordering operator P; however, 
fields at the same point do not commute and this fact must be taken into account. As a rule, one 
can make calculations in the abelian case; by rewriting the final result in a gauge invariant form, 
an expression also valid in the non abelian case is obtained. 
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SW q - q = -| (TrP^' f SS^izjF^zJ exp (ig f^dx? A^xty (2.28) 

where SS ttu (zi) = \{dz^8z" — dz^Szi) is the element of the surface spanned by z\{t). 
Then 

SW - r f { 

5(i\nW q7} )=i—^ = g 5S^{ Zl ){{F, u { Zl ))) (2.29) 

Wqq Jti 



and we may write 

*«^w» = K¥ (2 ' 30) 

(see App. A for a definition of ss ^/ z -, )■ The computation is similar for the case of Z2 with 
a minus sign of difference. 

Similarly it can be seen that 
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and therefore 
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From (1.13) we have then 
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-d a \d v D w {z x - z 2 ) - d p D up ( Zl - z 2 )}} . (2.34) 

Notice that obviously in the terminology of App. A we have C(z±, z^) = and then 

([(F»M)F P M))) - ((F,M))) ((F P M)})) SR = • (2.35) 
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By using the Coulomb gauge and the instantaneous approximation we have then 
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In a similar way we obtain also 
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(2.40) 



Let us now consider the confinement part of ilnWgg. From Eq. (2.23), taking into 
account that it™ 111 satisfies the appropriate Euler equation, one has (see App. B for details) 

fdu™ hl \ 
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where the subscript 1 indicates that the derivative is calculated in s = 1. A similar formula 
is valid for z 2 . Moreover as in the short range case 



<<J^(*i)i^(*i)>> " «i^(*i)» ((FM))) 



By using the straight line approximation one obtains 



LK 



0. 



(2.42) 



<7«*b*(*;)» 



LR 



a— + 0(v 2 
r 



a 1 



g ((F hk ( Zj ))) Lli = (rV - r k vT) + 0{v A ) 



rrij r 



(2.43) 
(2.44) 
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with j = 1,2 and 

/ \ LR ft 

g 2 [i(F hk (z 1 )F lm (z 2 ))) ~ (CM*i)» ((F lm (z2)))) = -*9js^) {{Fhk{Zl)))LR = 

= 0(v 2 ) . (2.45) 

Finally replacing (2.35)-(2.40) and (2.42)-(2.45) in (2.16) we obtain the terms in Eq. 
(1.3) involving explicitly spin. 

As concerns the Darwin-type terms one should evaluate g {{D v F vlx (x))) in Eq. (2.16). 
By using the results of this Section one can obtain by gauge invariance 



gf dz?((irF v ^z 1 ))) = f t t dt%&'((F VIM (z 1 ))) 

g£ dt (V((F„ » + ^((F^)))) = £ dtd h (- ^aj^ - a 



itj \ J Ju \ 6 r J r 

from which finally we have the Darwin terms as reported in Eq. (1.3). 

III. THREE-QUARK POTENTIAL 

Let us define the three-quark color singlet state (again for three quarks of different 
flavours) : 

^^^(yO^Cife)^^ (3.1) 



where the path-ordered gauge strings U are defined in Eq. (2.2), yu is defined as explained 
after Eq. (1.17), bi and di for i = 1,2,3 are the color indices, e^^ is the completely 
antisymmetric tensor in the color indices which puts the system of three quarks in a color 
singlet state and -4= is the normalization factor. 

The corresponding gauge invariant Green function can be written 

G(xi,X2,X 3 ,yi,y 2 ,y3) = ^a 1 a 2 a^b 1 b 2 b^ 

(0\TU asCs (x M ,x 3 )U a ^(x M ,x 2 )U a ^(x M ,x 1 )^ 3c3 (x 3 )^ 2c2 (x 2 )ip lci (x 1 ) 
^(yi)^^)^^)^ 161 ^!,^)^ 2 ' 2 ^,^)^ 363 ^,^)^) (3.2) 
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and we assume x\ = x 2 = x® = x° M = if, 2/1=2/2=2/3 = Vm = *ij t = U — U (t > 0). 
The integration over the Grassmann variables is again trivial and one can write 

1 / / \aibi 

G(xi, x 2 , x 3 , yi, y 2 , y 3 ; r) = ^e aia2a3 e blb2b3 ((U(x M , x 1 )Sf(x 1: y 1 \A)U(y 1: y M )J 

\a2b2 / \ a3&3\ 

U(x M ,x 2 )S 2 (x 2 , y 2 \A)U(y 2 , y M )) [U(xM,x 3 )S 3 (x 3 ,y 3 \A)U(y 3 ,y M )) ); (3.3) 



if some of the quarks are identical we have simply to sum over all permutations of the corre- 
sponding final variables. By performing even in this case the appropriate Foldy-Wouthuysen 
transformations and using (2.9) we find in place of (2.10) 

K(x 1: x 2 , x 3 , yi, y 2 , y 3 ; r\ x M , y M ) = 

/•zi(t()=xi rz 2 (t{)=X2 rz 3 (t ( )=x 3 



rz 1 (t i )=x 1 rZ2(t{)=X2 rz 3 [tf)=x 3 

= V[z uPl ] P[z 2 ,p 2 ]/ V[z 3 ,p 

Jzi(t i )=yi Jzo{U)=Y2 Jzi(ti)=v.i 

exp < i / dt y^ 



P? Pi ' 
p, • z,- — m,- - — I -z 

2rrij 8rrij 



-ee TT T s P exp J ig [_ dx» AJx) + — f dx^ 
3! " [ Jr- m^r, 



l3! ' 

(*$,(*) - ^^^ " ^ F ^) }) • ^ 

i^ being now the three-quark Pauli-type Green function. 

In (3.4) we have suppressed for convenience the color indices but have left trace of the 
tensors e^^a^b^b-i with the notation ee. As above T s denotes the chronological ordering 
for the spin matrices and P is the path-ordering prescription acting on the gauge matrices; 
Tj and Tj are defined as in Eq. (1.17) and following. Notice that the curve T made by the 
union of Ti, T 2 and T 3 is a closed three-branch loop which generalizes the Wilson loop of 
the two-body system. 

From now on one can proceed strictly as in Sec. II. We shall show that from (1.18) using 
(1.16) and (1.19) one can write 

3 

ee 1 
3! 



I 3 \ \ f tf 

-ee J| T s P exp . . . J ~ T s exp -i dtV 3q (z 1 ,z 2 ,z 3 ,p 1 ,p 2 ,p 3 ,S 1 ,S 2 ,S 3 



i=i 
and so the propagator K obeys the three-particle Schrodinger-like equation 
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(3.5) 



4k 

dt 



3 



p 



p 



. . v ,,? 2m.,- 8m^ 



y3g 



K 



(3.6) 



Again expanding the logarithm of the left-hand side, Eq.(3.5) turns to be equivalent to 
i\nW 3q + iJ2^j r dx»(s\ ((F ¥ (x))) - i^S\e lkr V ) ((F, r (x))) - 

7 = 1 J 3 J 



I , ,_ \ 1^ ig 2 



5m, 



((irF^x))))--^ 



2 r-f. vfinVfi^ 
n.i' J J 



:ij 



T s / dx» 



m,{x)FUx'))) - ({F ¥ (x))) ({F kff (x'))) 

with Wz q being defined as in Eq. (1.17) and now 

^{££{U j P[exp(igk i dx' i A lt (x))]}f[A]) 



(XOb O A tJj/ 



dtV 3q 



({f[A})) 



h(££{Uj Pexp(igfr dx^ A„(x))}) 



(3.7) 



(3i 



3. Contribution to the potential coming from i In W; 



■Ml 



Having in mind (1.19) let's evaluate 5 m i n in the same manner as we did for the two-body 
case. The quantity S'min is the area made by three sheet surfaces as described in Sec. I. Let us 
denote by Zm($) an arbitrary world line joining yu to xm and by Uj = Uj(s, t) the equation 
of an arbitrary sheet interpolating between the trajectories zf = Zj(t) and z\ f = z^(t); 
obviously Uj(0, t) = z^(t), Uj(l,t) = Zj{t) and u^(s,t) = t. 

Assuming that the minimum is taken in the choice of Uj (s, t) and of zjw(£) we can write 
in analogy with (2.23) 



S m - m = min V" / dt ds 

~{ Jti JO 



duj 



ds 



' duj 
~dt 



i 

2^ 2 



where now the index T, stands for transverse part of a vector with respect to 



ds 



d\ij 



ds 



Then, performing the straight line approximation (1.19) we have 



duf n 
ds 



Zj(t)-Z M (t) =Tj(t) 



(3.9) 



(3.10) 



(3.11) 
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duf 11 




dt 


and 


expanding in 


the velocities 






•Jmin 


-i: 


3 r 

dt Y, rj 

3=1 L 


with 


§ i = f j 









s±j(t) + (1 - s)z M (t) 



(3.12) 



1 ~ a&hi + ^MT, + ZjT, • Zjjfr.) + 



6 ^ 



(3.13) 



Taking into account this result and introducing (2. 18), (2. 20) and (2.21) in (1.18): 

2 a, 12 a, 



ilnW- 



3r/ 



jO te 



j<Z L "'J 



3^ + 23^^ + ^»^' 



+ 



+ aJ2' 



3=1 



1 ~~ 6^ jTj + * MTj + ijTj ' * MTj 



(3.14) 



where again r,-; = r~— r; = z,-— Zj. In the I configuration case the quantity % can be obtained 



deriving the equation Z)j =1 (r,-/r,-) — 0. We get J2j=iy-(3 — 



O -L f \ ft f\j fy, ft fy, f% \ ry !\> 



3 3 > 3 



n=^ hk - 



h~k\ -k 






j)z%j. Obviously in the II configuration case we have % = zi (M = quark/). 
Finally replacing Zj by Pj/rrij we obtain Eq. (1.8) 



4- Spin related potential 



As in the qq case we can write 



9((FM)) 



S(i\nW- 



■1<1) 



SS^(zj) 



(3.15) 



g 2 ((F^FpM))) ~ ((F^Zj))) ((F P M))) = W 



6 



-((F„„( Zj ))) (3.16) 



with j,i = 1, 2, 3 and adapt immediately the procedure used in the derivation of Eqs. (2.35)- 
(2.40) and (2.42)-(2.45) to the variation of each single quark world line separately. 
From the short range part of (1.18) we have then 

2 „ r r*i 



fj {(F t „,(zj)))™ = ^S*{J* dtii&D^Zj - zj - d^D up ( Zj - Zi )]zf+ 
+ / dtnilduD^Zj - z n ) - d^D up (zj - z n )]z^ 



(3.17) 



with j, i,n = cyclic permutation of 1, 2, 3 and 



g 2 ({(F, u { Zj )F pa {z t ))) ~ ((F^(zj))) m„{z t )))) 



SR 



;g 2 [d a d u D fip (z j - Zi) - d a d il D vp {zj - z t ) - 
- dpdyD^Zj - z^ + dpdpD^Zj - Zi)\ 



(3.18) 



for j ^ i; furthermore 



«iV(,,)i^(4)» " (CM%')» ((Fprty)) 



SR 



0. 



(3.19) 



By using the Coulomb gauge and the instantaneous approximation we have from Eqs. 
(3.17) and (3.18) 

2 



g{{Fok(zj))) 



fy-4 n> fy, Fb * 

u V ' ji jn / 



SR 
'3)11 q" s I „3 '3 I ' 



(3.20) 



O II li I jj^ 

3 m r • -'" J nFn ' 



"n ' Jn 



(j, i,n — cyclic permutation of 1, 2, 3) and 



g 2 ((F hk ( Zj )F lm (z t )}} - ((F hk ( Zj )}} {(F lm ( Zi ))) 



SR 



-1^-"){W 






°i u h 



(3.21) 



/ rfcm , -fe -to 



-^* } 



j - ' 



S hl + f h r l - 

3 l J' 1 






J'' 



S kl + r k r l - 



(3.22) 



where j ^ i- 

Let us now consider the long range part of Eq. (1.18). Notice that, due to its definition, 
an infinitesimal variation of the intermediate point world line z^ = zj^(t) leaves the quantity 
% In Wz q unchanged. Then, in the evaluation of the functional derivatives (3.15) and (3.16), 
we can treat such world line as fixed. From Eq. (3.9), taking into account that w™ 11 satisfies 
the appropriate Euler equation, one has 
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9({F,u(zj))) 



LR 



a 






Zj v 



X < — Z 



I du mm s 

2 / ^"j 



9s 



J / 1 



'du™ n \ 



'duf m \ 
, 9s j )i. 



x 



where the subscript 1 indicates that the derivative is calculated in s = 1. 
Moreover 



By using the straight line approximation one obtains 

9 «*b*te)» LR = a 



LK 



0. 



r j 



(3.23) 



(3.24) 



(3.25) 



with j = 1,2,3 and 



g ((F hk ( Zj ))) LR = — -(rfrj - r*pf) 



LR 



(3.26) 



(? 2 ((F ftfe (^)^m(^)» " (CM^)» ((Fim(zi))) = ig 



SS lm ( Zi ) 
0(v 2 ) 



(CM^)» 



LR 



(3.27) 



with j j^ i. 

As concerns the Darwin-type terms we must evaluate g {{D^' v F v ^{zj))) for j = 1,2,3 
with D^ u = d^ u + igA u . This can be done by using the same line of derivation as for the 
two-body case. 

At this point we may derive the spin dependent potential: after some calculations using 
Eqs. (3.19)-(3.22) and (3.24)-(3.27) in Eq. (3.7) we get Eq. (1.7). 

IV. ADDITIONAL CONSIDERATIONS 

To complete the work some additional considerations and remarks are necessary. 
1) In Sec. Ill we have assumed different flavours for the three quarks. If two or three 
quarks are identical we must identify the corresponding operators in Eq. (3.1) and add the 
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appropriate normalization factor -4= or -4=. Correspondingly we have to replace the right- 
hand side in (3.3) with the corresponding sum over the permutations of the final identical 
particles divided by the combinatorial factor (2! or 3!). For what concerns the potential this 
means only that we have to equate the masses of such particles in (1.7) and (1.8). Similarly 
if in Sec. II the quark and the antiquark have the same flavour, a new term must be added 
to the last member of (2.1), which is obtained by interchanging the roles of y\ and x-i. Such 
"annihilation" term is not properly of potential type but can be treated perturbatively [10]. 
2) Concerning the ordering problem in (1.4) and (1.8) we recall that there are two 
independent possible prescriptions for a quantity quadratic in the momenta: 
the Weyl prescription 

{f(r)p h p k U = \{{m,p h },p k } (4.1) 



and the symmetric one 



{/(r)pV}s = ^{/(r),pV}. (4.2) 



As well known, in the path integral formalism the ordering prescription corresponds to the 
specific discretization rules used in the definition of it. 

In Ref. [5] we adopted a somewhat ad hoc rule corresponding to the ordering { } or d = 
|{ }w + §{ }s- Such rule was motivated by the fact that it enables a by part integration at 
the discrete level which was necessary to eliminate a dependence of the potential from the 
acceleration. 

Notice however that the limit procedure used for the definition of the path integrals in 
(2.9) or (2.10) is not at our choice but it is a consequence of the corresponding procedure 
used in the definition of the field functional integration in (2.1). To see what is the correct 
prescription, let us assume a definite lattice with spacing e in the time direction and a 
spacing a in the space directions. Let us consider the corresponding discrete counterpart of 
(2.1), written according to the usual rules for gauge theories [11] and perform the integration 
of the fermionic fields at this discrete level. Then in place of (2.10) we arrive at an equation 
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in which not only the time integral is replaced by a sum over the appropriate discrete times 
t s , but for every t s even the integrals on Zi and z 2 are replaced by the sum over all the sites 
of the lattice corresponding to that time coordinate. Finally Eq. (1.12) has to be interpreted 

as 



Wqq — 



j U dU n'ne tS[u] TrP I] U *r (4.3) 

{n',n} {r,r'}er 

where U n > n denotes the element of the color group associated to the link between the con- 
tiguous sites n and n' and the product is extended to all such links or to all links laying on 
the curve T. Since in turn U n > n can be interpreted as exp ig(x n > — x n )^A^( Xn ' 2 X " ) , we see 
that, after having explicitly integrated over U (and so used (1.13), (1.14)) and performed 
the limit a — » 0, we are left with the discrete form of an ordinary path integral with 

X hk (r) P ^ — > X hk ( r -^^) P t P % . (4.4) 

Eq. (4.4) does correspond to the Weyl ordering, as indicated in (1.4) [12]. Notice, how- 
ever, that a different ordering in (1.4) would bring simply to additional Darwin-like terms 
which in practice can be nearly completely compensated by a readjustment of the potential 
parameters. Similar arguments apply to Eq. (1.8) in the 3g case. 

3) To clarify the connection between the qq potential and the relativistic flux tube model 
[2] it is convenient to neglect the spin dependent terms in (2.10) and replace the 1/m 2 
expansion by its exact relativistic expression. We have 



i^(xi,x 2 ;yi,y 2 ;tf-ti) 



/ Z>[zi, p 2 ] / V[z 2 , p 2 ] exp < i j dtJ2(Pj ■ Zj - \fm) + pf) 



>u j= i 



+ \nW qq -\ . (4.5) 



Further by taking advantage of (2.22) and of the first step in (2.27), after expanding again 



the exponential in (4.5) around the values Pj = mZj/Jl — z| and performing the integration 
in the gaussian approximation (semiclassical approximation), we can write (see however in 
this connection ref[12]) 
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Zf(xi,x 2 ;yi,y2,tf-ti) = 

/ P[zi]A[zi] / £>[z 2 ]A[z 2 ]exp ji / dtL(z 1 ,z 2 ,z 1 ,z 2 



(4.6) 



with 



3=1 6 T 

— ar ds[l — (szit + (1 
Jo 



l--(5 hk +f h f k )z^z. 



+ 



S)Z 2T 



,211/2 



(4.7) 



In (4.6) T>[z] denotes the usual nonrelativistic configurational measure and A[z] is a deter- 
minantal factor which has to be be considered part of the relativistic measure. Formally we 
can write 



V[z] 



(e = time lattice spacing) and 



m \ 3 / 2 t-t 
2wk) V 



1 



2-Kit) y ' 



(U 



A[z] = mdet -(l-z 2 ) 1/2 (S hk -z h z k ) 
\ \_m 



(4.9) 



What we want to stress is that Eq. (4.7) is identical to the center of mass lagrangian for 
the relativistic flux tube model. This is consistent with what already observed at the 1/m 2 
order in Ref. [2]. 

4) In phenomenological analysis the following long range static potential of two-body 
type has been often adopted [7] for the 3q case: 



V. 



^t = ^(n 2 + r 23 + r 31 ] 



(4.10) 



with a corresponding spin dependent potential again of the form (1.10) 



K LR 






(4.11) 



The factor 1/2 in Eq. (4.10) is motivated by the fact that when two quarks collapse they 
become equivalent from the colour point of view to an antiquark and V st l t: V S( ± must reduce 
to K-tat and V^. It has been shown that (4.10) and (4.11) produce a spectrum very close 
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to that obtained from (1.6) and (1.10). From a numerical point of view the use of (1.6) 
amounts to replace the factor 1/2 in (4.10) by a factor of the order of 0.54-0.55 [7,8]. Notice 
however that from a fundamental point of view (4.10) has no clear basis. 

5) As we already mentioned in lattice gauge theory the area law is obtained under 
the approximation in which the quantity Mf(A) is replaced by 1, the so called quenched 
approximation. Our entire treatment has to be understood in this perspective. Then the 
effect of virtual quark-antiquark creation should be introduced as a correction at a later 
stage. Various attempts have been done in this direction: see references [13,14]. 

6) In Eqs. (1.6)-(1.8) V s ^ has been written in the coordinates ri,r-2,r3 and ri2,r23,r3i 
which provides a symmetric form of the potential. For numerical computations it may be 
useful to express such variables in terms of a system of independent coordinates like the 
Jacobi coordinates (see [3] and [8]). 

V. CONCLUSIONS 

In conclusion we have strongly simplified the derivation of the quark-antiquark potential 
as given in [4,5]. We have shown that, once the assumptions (1.13) and (1.14) are done, the 
existence of a potential follows when one performs the instantaneous approximation (1.16) 
and the straight line approximation (1.15). We have also corrected the ordering prescription. 

Due to the above simplifications, the method has been extended without difficulty to 
the three-quark case, where the relevant assumptions and approximations are (1.18), (1.16) 
and (1.19). As a result a 3q spin dependent potential has been consistently obtained in 
the Wilson loop context. This coincides with the one already proposed by Ford under 
an assumption of scalar confinement. It has also been evaluated the 0{l/m 2 ) 3q velocity 
dependent potential which is new at our knowledge. 

Notice that in both the qq and 3g cases the spin independent relativistic corrections ob- 
tained by us differ from those resulting from the mentioned assumption of scalar confinement 
and seems to agree better with the data [14,10] (for the difficulties of the scalar confinement 
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hypothesis see also [15,16]). 

Finally we have seen that if in Eq. (2.10) we replace the kinetic terms by the exact 
relativistic expressions and neglect the spin dependent terms, if further we use the first step 
in Eq. (2.27) (without any velocity expansion) and perform the momentum integration in 
the semiclassical approximation, we obtain the lagrangian for the relativistic flux tube model 
[2]- 

APPENDIX A: 

Concerning the definition of the symbol ssi t, z \ , let us consider a functional of the world 
line 7 of the form 

F[ 7 ] = A+ f dz»B,{z) + l - I dz» f dz'P C pp (z, z') + ... (Al) 

with C pp (z, z') = C pp (z, z') etc. Let us denote by z = z(X) the parametric equation of 7 and 
consider an infinitesimal variation of such curve z(X) — ► z(X) + Sz(X). We have 

SF= f (sdz^B^z) + dz» dB ^ 5z u \ + 



l dCfj f {z 1 z > ) 

./- V""" ' '"""' " ' dz v 

'dBJz) dBJz) 



+ J dz'o J Udz^c^z, z r ) + dz ^'":}:r' s-j- \ - - . . . 

= - f(dz"5z' / -dz u 5z^ 

f jjpf 90 ^^) dC vp {z,z') \ 
+ i 7 I dz" dz» + --- 



dz u dz» 



(A2) 



where a partial integration has been performed at the second step. Assuming 5z(X) different 
from zero only in a small neighbourhood of a specific value A of A we write 



SF dBJz) c)B v $)_ + r dz , p f dC w (z,z') _ dC vp {z,z'Y 



+ ... (A3) 



SS^(z) dz v dz» }-< \ dz v dz» 

with ~z = z(X). Furthermore if we consider a second variation Sz(X) different from zero only 
in a small neighbourhood of a second value A ^ Awe have 
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5 



5F 



SS^(z) 



(dz' p 5z' a - dz' a 5z' p ) 



d fdCUz,zf) dC up (z,z>) 



dz 1 ' 



dz v 



dzP 



d dC^(z,z') 



dz'p 



dz h 



dC ua (z,z') 
dzP 



+ 



and 



5 2 F 



d 2 C, p (z,z') d 2 C up (z,z') 



5Sp°(z')8S^(z) 



dz' a dz u 



dz' u dzP 



d 2 C pa (z,z') | d 2 C V(J {z^) | 



(A4) 



(A5) 



dz'Pdz u dz'PdzP 

Notice that the assumption A' 7^ A is essential to make the definitions non ambiguous. The 
case in which a 5(z — z') term occurres in C(z, z') must be treated as a limit one. In practice 
this amounts to say that (A3) and (A5) hold true even in this case. 



APPENDIX B: 



In applying (A3) and (A5) to i\nW qq - [17] and ilnW^q and specifically to Eqs. (2.23) 
and (3.9) it is convenient to think of t as an arbitrary parameter on the same foot of A. 
It is obviously understood that at the end t is identified with the ordinary time by setting 



z°At) = t and u°(s,t) = t (this in the qq case). In particular rewriting Eq. (2.23) as 



ilnW, 



LR 



a 



<7 ■ 



rr I dt dsS(u r 

hi Jo 



(Bl) 



with 



S(u) 



' du p du,, \ ( dup du. 



' du p du 



dt dt J \ ds ds J \ dt ds 

one can notice that the equation of the minimal surface u p = a 

Euler equations 

d dS d dS 





(B2) 
' ' s, t) is the solution of the 



ds dm *d (%) 



(B3) 



satisfying the contour conditions u M (l,i) = Zi(t), u p (0,t) = z^t). Then considering an 
infinitesimal variation of the world line of the quark 1, z p {t) — ► Zi(t)+Szi(t), even u^ lin (s, t) 



must change, Mm in (s,t) 



u,. 



(s, t) + 5uP(s, t) and one has 
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5(i\nWJf) = a J' dt f ds 



OS d . u dS d s u 



d (*g)d8 



d(^ 



dt 



rtf 

a dt 



dS 



d( d i£ 



Su 1 " 



(B4) 



s=l 



where Szi(t) was assumed to vanish out of a small neighbourhood of a specific value of t. 
Finally taking into account that 



8u»(l,t) =Sz?(t) 



dt 



%(t) 



one obtains 



5(t\nW^)=ajytj^— 



■ 2 fdu?" 

h { ds 



+ 



1 



ds 



Z\Z\ V 



1 



1 /•** 

—a \ dt (dziSzi — dz\5z 

2 Jt; 



and so Eq. (2.41). 



" x ~0 



'du™ m \ . (du 

Z\ u — 



ds 



i 



x < -z. 



'du 1 - 



ds 



+ 



Zl 



ds 



'du 1 



bz\ 

Zlfi 



ds 



2\ ~2 



(B5) 



(B6) 
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FIGURES 
FIG. 1. The two types of configurations of the three-quark system. 

FIG. 2. Generalized Wilson loop for the quark-antiquark system. 

FIG. 3. The analogous of the Wilson loop for the three-quark system. 
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